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Abstract 



X : When you measure an observable, A, in Quantum Mechanics, the state of the system changes. 

This, in turn, affects the quantum-mechanical uncertainty in some non-commuting observ- 
able, B. The standard Uncertainty Relation puts a lower bound on the uncertainty of B 
in the initial state. What is relevant for a subsequent measurement of B, however, is the 
uncertainty of B in the post-measurement state. We make some remarks on the latter prob- 
lem, both in the case where A has a pure point spectrum and in the case where A has a 
continuous spectrum. 
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1. Introduction 



The Uncertainty Principle is one of the signature features of Quantum Mechanics. In popular 
accounts, it is often described as the principle that measuring some observable "disturbs" 
the system and that this has implications for subsequent observations. As we shall review 
in £j2j the Uncertainty Principle does not say anything about successive measurements. In 
fact, formulating a precise statement about the uncertainties in the result of successive 
measurements (of non-commuting observables) has received, perhaps, less attention than it 
deserves. 

The character of such a statement depends a great deal on whether the observable(s) 
in question have discrete or continuous spectra. We start with the former case, in ^2], and 
tackle the latter in §H When the spectrum is continuous, we need to include a finite detector 
resolution as part of what it means to measure that observable. This has dramatic (and, to 
our knowledge, quite novel) implications for the results of successive measurements. In §4.11 
we find, for instance, that if one first measures x (with detector resolution, a x ) and then 
measure p (with detector resolution o~ p ), the product of the uncertainties of the measured 
values is bounded below by 

( Ax ) measured ( A P) measured > \ (} +\J 1 + 4c ^p) 

which contrasts quite strikingly with the lower bound on the product of quantum-mechanical 
uncertainties in the initial state: (Ax) 2 p (Ap) 2 p > |. In §4.2[ we extend this formalism to the 
notion of "joint measurements" of non-commuting observables. 

Our interest in these questions was sparked by an attempt to understand pQ, which is one 
attempt to study this problem of successive measurements. In §3 we make some remarks 
about the relation between our two approaches. In §5.11 we derive an "intrinsic" expression 
for the "disturbance" rj(B), of some observable B, which results from measuring some other 
observable, A. In §5.21 we turn to Lund and Wiseman's refinement of Ozawa's proposal, 
involving a "weak measurement" of A. We show how this can be cast in terms of a family 
of Positive Operator- valued measures (POVMs), and find the corresponding expressions for 
rj(B) and e(A). 

Our interest in [1] was, in turn, inspired by a recent experimental attempt [2] to test 
these ideas. In the end, our results shed little light on the latter experiment. But we believe 
them to be of sufficient interest to stand on their own. 

While this paper was nearing completion, the preprint [3] appeared, which discusses 
related issues. 

The previous version of this manuscript used a less "efficient" POVM for measurements 
of x and p (essentially, the POVM for joint measurements, of §4.21 where one discards 
the result for one of the jointly-measured observables). The inefficiency of that choice was 
pointed out to us, in private communication, by A. Di Lorenzo, whose paper |3] contains 
results overlapping with those of §U 
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2. Measuring an Observable with a Pure Point Spectrum 
2.1. General setup 

A quantum system is described in terms of a density matrix, p : % — > H, which is a 
self-adjoint, positive-semidefinite trace-class operator, satisfying 

Tr{p) = 1, Tr(p 2 ) < 1 
In the Schrodinger picture (which we will use), it evolves unitarily in time 

p(t 2 ) = t/(t 2 ,ti)p(ti)f/(t 2 ,ti) _1 (1) 

except when a measurement is made. 

Consider a self-adjoint operator A (an "observable"). We will assume that A has a pure 
point spectrum, and let be the projection onto the i th eigenspace of A. 

When we measure A, quantum mechanics computes for us 

1. A classical probability distribution for the values on the readout panel of the measuring 
apparatus. The moments of this probability distribution are computed by taking traces. 
The 77 th moment is 

(A n ) = Tr(A n p) 

In particular, the variance is 

(AA) 2 = Tr(A 2 p) - {Tr(Ap)f 

2. A change (which, under the assumptions stated, can be approximated as occurring 
instantaneously) in the density matrix [5], 

p^=p(p,A) = J2Pi iA) pP l iA) (2) 

i 

Thereafter, the system, described by the new density matrix, p, again evolves unitarily, 
according to ([T]). 

The new density matrix, p, after the measurement , can be completely characterized by two 
properties 

frequently, one wants to ask questions about conditional probabilites: "Given that a measurement of A 
yields the value Ai , what is the probability distribution for a subsequent measurement of ..." . To answer 
such questions, one typically works with a new ("projected") density matrix, p' = ^P[ A ^ pP[ A \ where the 

normalization factor Z = Tr(p[ A ^ p) is required to make Tr(p') = 1. The formalism in the main text of this 
post is geared, instead, to computing joint probability distributions. 
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1. All of the moments of A are the sam^l as before 

(A n ) = Tr(A n p) = Tr{A n p) 

(In particular, AA is unchanged.) Moreover, for any observable, C, which commutes 
with A ( [C, A] — ), its moments are also unchanged 

(C n ) = Tr{C n p) = Tr{C n p) 

2. However, the measurement has destroyed all interference between the different eigenspaces 
of A 

Tr([A,B]p) = 0, VB 

Note that it is really important that we have assumed a pure point spectrum. If A has a 
continuous spectrum, then you have to deal with complications both physical and mathe- 
matical. Mathematically, you need to deal with the complications of the Spectral Theorem; 
physically, you have to put in finite detector resolutions, in order to make proper sense of 
what a "measurement" does. Later, we will explain how to deal with those complications. 

Now consider two such observables, A and B. The Uncertainty Principle gives a lower 
bound on the product 

(AA) p (AB) p >l\Tr([A,B]p)\ (3) 

in any state, p. 

The proof, generalizing the usual proof presented for "pure states", is as follows. Let 
S=(A- (A) p l) + e ie a (b - (B) p l) for a G R. Consider 

Q(a) = Tr{SS j p) 

This is a quadratic expression in a, which is positive-semidefinite, Q{a) > 0. Thus, the 
discriminant must be negative-semidefinite, D < 0. For 9 = ir/2, this yields the conventional 
uncertainty relation, 

{AAf p {ABf p >- A {TrmB]p)f 

For 9 = 0, it yields 

{AAf p {ABf p > (^Tr({A, B}p) - Tr{Ap)Tr{Bp)^ ' (4) 

which is an expression you sometimes see, in the higher- quality textbooks. 

2 More precisely, the probability measure on Spec(A) is unchanged (see M3. 1[) . This is an important 
distinction when A is an unbounded operator, as the moments of A generically don't all exist, and so can't 
be used to characterize the probability distribution. 
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As stated, ([3]) is not a statement about the uncertainties in any actual sequence of 
measurements. After all, once you measure A, in state p, the density matrix changes, 
according to (j2j), to 

P(P,A) = J2P 1 {A) PP 1 (A) (5) 

i 

so a subsequent measurement of B is made in a different state from the initial one. 

If we are interested in the uncertainties associated to successive measurements, the ob- 
vious next thing to try is to note that, since the uncertainty of A in the state A) is the 
same as in the state p, and since we are measuring B in the state p, we can apply the Uncer- 
tainty Relation, (j3J) in the state p, instead of in the state, p. Unfortunately, Tr([A, B]p) = 0, 
so this leads to an uninteresting lower bound on the product the uncertainties 

(AA)(AB) = (AA) p (AB) p >0 (6) 

for a measurement of A immediately followed by a measurement of B. 
To get a (slightly) more interesting lower bound, let 

B = J2 P i A)BP i A) (7) 



and 



We easily compute 



and hence 



Mi = P„ (A) B (l- P (A) 



(AB)|= (AB) 2 p + J2 Tr ( M * M iP) 



(AA) 2 p (AB) 2 p = (AA) 2 p (AB) 2 p + (AA) 2 p £ Tt^mIp) 

i 

Since [B, A] = 0, the first term is bounded below (see (J3|)) by 



(AA) 2 p (AB) 2 p > (\Tr{{A, B}p) - Tr{Ap)Tr{Bp)\ (9) 



The second term is also positive-semidefinite. So, the product of the measured uncertainties, 
for successive measurements of A then B, obeys 



(AAf(AB) 2 = (AA) 2 (AB)l 



> 



-Tr({A, B}p) - Tr(Ap)Tr(Bp) J + (AA) 2 p £ TriM^p) 



(10) 



4 



It is conventional to define a quantity 



Tr 



((£? — B) + 3-B 2 )p 



called "the disturbance in B" (due to the measurement of A). Here B is defined as in (J7J) 
(and similarly for B 2 ), and the trace is taken in the initial (pre- measurement) state, p. The 
combination B 2 — B 2 , which appears in (fTTT) . is exactly the positive semi-definite operator, 
J^^MjM/, which appeared in (fTUl) . 



2.2. The two-state system 

For the classic case of a 2-state system, with A = J x and B = J y , we see that B = 0, and 
the product of uncertainties is entirely given by the second term of ([8]). 

The most general density matrix for the 2-state system is parametrized by the unit 3-balJH 



P 



1 



1 + a ■ a) 



a ■ a < 1 



The points on the boundary S 2 = {a ■ a = 1} correspond to pure states. 

Upon measuring A = J x = ^a x , the density matrix after the measurement is 



P 



and, for a subsequent measurement of J, 



- (1 + a x a x ) 



(AJ x ) 2 p (AJ y )l 



16 



A- < 



as "predicted" by ©. 

More generally, consider B = b ■ J. The quantum-mechanical uncertainty in the initial 
state is given by 



(ABY 



1 



(a- by 



'P 4 

After measuring A = J x , the quantum-mechanical uncertainty in the state p is given by 



{AB%=^[b 2 -a 2 x b 2 x ] 



so the change in the uncertainty, 



(AB%-(AB) 2 p = - 



(a-by 



2l2 



(12) 



can be arbitrarily negative. If we fix b 2 = 1, and take a = t=(1, 1, 0), corresponding to the 
pure state 



3 More generally, for an n-dimensional Hilbert space, "H, the set of density matrices is an (n 2 — 1)- 
dimcnsional convex set S C R"" _1 . The (n 2 — 2)-dimensional boundary of S consists of density matrices on 
proper subspaces T~L' C T~L. The set of pure states is a closed (2n — 2)-dimensional subset of the boundary 
consisting of density matrices of rank-1. This 11-1 consists of the extreme points of S [BJ. 
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p 



1 

givr/4 



e -»7T/4 
1 



then the minimum for (fT2|) is achieved for 6 = J=(l, — 1, 0) (i.e., £> = -^(-^ — 



(AB)l-(AB) 2 
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(13) 



Nevertheless, the disturbance, 



1 



> 



Note however that, for the 2-state system, r)(B) depends only on the relative orientation of 
A = J x and B = b ■ J and not at all on the state, p, in which it is measured. We will return 
to this point in §5.11 

Before we turn to the measurement of observables with continuous spectra, we need to 
make a disgression on the Spectral Theorem, Projection- Values Measures, and the more 
general concept of Positive Operator- Valued Measures. 



3. Project ion- valued and Positive operator-valued measures 
3.1. PVMs and observables 

For an observable, A, with spectrum X = Spec(A), it is rather awkward (and, in the 
unbounded case, generically impossible) to describe a probability distribution implicitly, in 
terms of its moments, (A n ) . Instead, we would like to directly give a formula for the desired 
probability measure on X. 

That is what the formalism of projection- valued measures (PVMs) affords us (see, e.g., 
[7]). A projection- valued measure on X is a rule which assigns to every Borel subset S <Z X 
a projection operator tts on H satisfying the obvious properties 

• 71 X = 1. 

• 7T0 = 0. 



• TTSjUSau... = 7r Sl + 7T 52 + . . . whenever Si n Sj = 0. 

The Spectral Theorem is the statement that there is a 1-1 correspondence between self- 
adjoint operators, A, with Spec(A) = X and projection- valued measures on X. 
For % = £ 2 (r) and A = x, the corresponding projection- valued measure is 

(tt s ^){x) = \ . (14) 

I (J otherwise 
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For A = p, the corresponding projection-valued measure (slightly schematically) is 



jdy J ^e*<*-»ty(v) 



With the Spectral Theorem in hand, it is easy to give a more explicit description of the 
probability measure on X = Spec(A). It is the measure which assigns to each Borel subset, 
S C Spec(A), the probability 



Unfortunately, when X is a continuous spectrum, there's no generalization of von Neu- 
mann's formula ([2]) for the change in the state, due to the measurement. Indeed, no formula 
is possible until we introduce the notion of a finite detector resolution. But, in order to 
do that, we need a more expansive notion of what is "measurable" in Quantum Mechanics, 
which brings us to the concept of a Positive Operator- Valued Measure. 

3.2. POVMs and "noisy" measurements 

The Spectral Theorem establishes a correspondence between self-adjoint operators and projection- 
valued measures. Given a state, p, the projection- valued measure gives us a probability 
measure on X = Spec(A). But a projection- valued measure is not the most general way to 
manufacture such a probability measure on X, from a state, p. 

In discusssing observables with a pure point spectrum, we have assumed that our mea- 
suring apparati are "ideal." That is, we have assumed that the measured value of A was 
an accurate reflection of the quantum state, p, of the system. One extension of this anal- 
ysis would include a discussion of imperfect measuring apparati. For instance, even with 
a discrete spectrum, our measuring apparatus might be unable to distinguish between two 
closely-spaced eigenvalues Aj, Xj (this is a particularly pressing concern when the spectrum 
of eigenvalues has an accumulation point). 

In this case, the measurement need not destroy the interference between these eigenspaces. 
To account for this, we replace the pair of projections P^ A \Pj A ^ with a single projection 

Pif^ = + Pj A ^ i n von Neumann's formula for the change in the state as a result of 
the measurement. 

More insidiously, the measuring apparatus may not be perfectly classical in its behaviour. 
The different values for the pointer variable(s) may not decohere. In that case, the "mea- 
surement" produces an entangled state of the measuring apparatus and the system under 
study. For some purposes, this is actually an advantage and much current effort is devoted 
to exploiting the possibilities this affords. 

To account for these and other, more general, effects, one must leave the world where 
observables are self-adjoint operators (equivalently, projection-valued measures) and enter- 
tain a wider class of observables (see, e.g., [8] for an introduction), corresponding to positive 
operator- valued measures (POVMs). 

Neumark's Theorem guarantees that a POVM can always be represented as a PVM on 
some larger Hilbert space. When the resulting PVM has a pure point spectrum, the POVM 
is a countable collection of self-adjoint positive-semidefinite operators, {i^}, satisfying 



P(S) = Tr(TT S p) 



(15) 



7 



E^ 1 

i 

(Here Spec(Fi) C [0, 1], whereas projection operators had Speciji) = {0, 1}.) The probabil- 
ity of measuring the value Aj is given, in precise analogy with ( Fl5l) . by 

P(A<) = Tr(F,p) 

Unfortunately, even in this simplest case, the replacement for (J2j) is not determined by 
the POVM alone. We can write 

Fi = L\U (16) 

and (|2J) is replaced by 

Paftcr = P = ^ L iP L l ( 17 ) 

i 

But there are an infinite number of solutions to (flEl) . and hence an ambiguity in (117p . Most 
authors take the extra data of a choice of Lj, satisfying fll6p . as part of the prescription of 
the POVM, and we will follow that tradition. 

We will find that measuring — with finite detector resolution — an operator with a 
continuous spectrum, can be cast as a particular class of POVMs. 

4. Measuring Unbounded Operators with Continuous Spectra 
4.1. Successive measurements of x and p 

Let's go straight to the worst-case, of an unbounded operator, with Spec(A) = R. Such an 
operator has no eigenvectors at all. What happens when we measure such an observable? 
Clearly, the two conditions which characterized the change in the density matrix, in the case 
of a pure point spectrum, 

1. The probability distribution for measured values of A has moments given by 
Tr(A n p) = Tr(A n p). 

2. Tr([A,B]p) = 0, Vfi 

are going to have to be modified. The second condition clearly can't hold for all choice of 
B, in the continuous-spectrum case (think A = x and B = p). As to the first condition, the 
probability distribution, for the measured values of A, depend on the detector resolution. 
But the state p knows nothing about that. Instead, we will argue that the probability 
distribution for the measured values of A is given by a particular POVM. 

To keep things simple, let's specialize to % — £ 2 (R) and A = x. Imagine a detector 
which can measure the particle's position with a resolution, o~ x . We expect this detector 
resolution to add, in quadrature, to the quantum-mechanical uncertainty, of x in the state 
p, to give the uncertainty of the measured value of x. 
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A POVM which implements this is the following. Define the positive semi-definite oper- 
ator 

Fx = Ll L X0 (18) 

where L XQ is the operator 

L X0 : ip(x) H> f(x - XQ)i>{x) (19) 
If the function / satisfies the normalization condition 



1 = J du\f(u)\ 2 (20) 
then the F xo obey the completeness relation 

J dx F X0 = 1 (21) 
and define a POVM which assigns to a Borel subset S C Spec(x) the positive operator 

S m- JF S = [ dx F X0 (22) 



Because it simplifies several of the formulae which follow, we will always assume that this 
"acceptance function" is an even function: f(u) = f(—u). 

When we measure x, with this finite- resolution detector, the density matrix changes, 
according to the generalization of von Neumann's formula, (fTTj) : 

P 



J dx L XoP Ll Q (23) 



If p is represented by the integral kernel, K(x,y), then the state after the measurement, p, 
is represented by the kernel 

K(x, y) = J dx f(x - x )J(y - x ) K(x, y) (24) 

If the acceptance function, f(u), is sharply-peaked near u — 0, then f )24^) has the desired 
effect of suppressing the off-diagonal elements of K(x,y). 

The quantum- mechanical probability distribution for x is unaltered by (|23|) . Indeed, all 
the moments 

Tr(x fc p) = Tr(x fc p) 

and, in particular, (Ax)| = (Ax) 2 p . However, the probability distribution for the measured 
values of x is given, not by the PVM (114)) . but by the POVM f |22|) . which is sensitive to the 
detector resolution 

2 



( A %) measured = J dx Q X 2 Q Tr(F xo p) - dx X Tr (F xoP ) 



(25) 



< = J duu*\f{u)\ 2 (26) 
Thus we have the two properties we sought: 



where 

2 
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(|24|) expresses our expectation that measuring x suppresses the off-diagonal elements 
of the density matrix, p, after the measurement: K(x, y) — > for \x — y\ 3> a x . 



• ( 125|) expresses our expectation that the detector resolution, a x should add in quadrature 
with the quantum-mechanical uncertainty to give the measured uncertainty. 

While quantum-mechanical probability distribution for x was unaffected by the measure- 
ment, the probability distribution for p is altered by ( 1231) . The first few moments are 

Tripp) = Tr(pp) 
Tr(p 2 p) =Tr(p 2 p)+ri(p) 2 1 ' 

where 

V (p) 2 = - [ duf(u)f"(u)= [ du\f'(u)\ 2 



In particular, the quantum-mechanical uncertainty increases as a result of the measurement 

(Ap)| = (Ap)5 + V (p) 2 (28) 
An example, which will prove to be most useful, is a Gaussian acceptance function 

f(u) = ^— e -» 2 l^ (29) 

K(x, y) = e- {x - y)2/8cT *K(x, y) (30) 
v(p? = ± (3D 



For ( 1241) . we obtain 
and 



The advantage of the Gaussian acceptance function is that, for given detector resolution, a x , 
it minimizes the change ( 1281) in the quantum-mechanical uncertainty in p. 

Of course, the quantum-mechanical uncertainty in p is not quite what we measure. We 
need to impose a finite detector resolution for p, as well. Again, this corresponds to a POVM. 
Let us take a "Gaussian" detectoi0 

L k0 : H. [j=j=\ j dye-^-^ 2+ ^-y^(y) (32) 

and 

F fc0 = 4 o L fco : ^(x) ^ | dy e-^C-v) 8 -^^) ^) (33) 
which obey the completeness relation 

dkn ~ . 



4 This is just the momentum space version of the Gaussian detector we described for measuring x. 
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The uncertainty in the measured value 



(Admeasured = J ^ k l Tr G^oP) ~ (J ^ k Tr(F ko p) 



2 



(35) 



= (Ap)f+a 2 

Using (EHD, 



(Ap)Lasurcd = (Ap)|+^ 



V 

2 , „/„\2 , 2 



where, in the last line, we used that r](p) 2 is minimized for the Gaussian detector, (l3Tj) . 

From this, we see the product of uncertainties — for successive measurements, first of x 
then of p — obeys 

(Ax)Lasurcd(Ap)Lasu r ed = ((Az)J + <7*) ((A P )J + 7?( P ) 2 + 

> \ + (^2 + °j) ( A *)p + ^( A P)p + °M (37) 

where we used that (Ax) 2 (Ap) 2 > |. The last inequality is saturated by choosing p to be a 
pure state, consisting of a Gaussian wave packet, with very carefully chosen width: 

One can consider other acceptance functions in (j!9p . One nice choice is 

f(\_ f 2 _ ^ 

1+ cosh(^)' " 2coth 2 (a6)(a6coth(a6) - 1) 

cosn(ao) v / v \ / / 

which asymptotes to a square pulse (supported oniG [—b, b] ) in the limit a —> oo. Doing 
the requisite integrals is a little harder, in this case, but still eminently doable. We easily 
compute 

V n {x,y) = G{x-y)(-i^-\ 5(x - y) 



dx 



where 

G(u) = I dwf(w)f(w — u) 



asinh(a&) (b - u/2) sinh a(b + u/2) - (b + u/2) sinha(6 - u/2) 
2(ab cot\i(ab) — 1) sinh(mt/2) sinha(6 + u/2) sinha(6 — u/2) 

1 ° 2 ( 1 3 ^ u 2 + 0( 4 ) 

~ 12 Va6coth(afe) - 1 ~ sinh 2 (a6) / lU j 
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and hence 



P = P, P 2 = P 2 + 7/0) 2 l, r](p) 2 



6 \ab coth(ab) — 1 sinh (ab) 

Similarly, 



X = X, x 2 = x 2 + cr 2 1 



where 

2 [ . 2// n2 ,2/1 2a6cosh(a6) vr 2 



a 2 = / duu z f{uY = b l 1 - - - r/ir ^ in \ + 



3(a6cosh(a6) — sinh(afe)) 3(ab) 2 

These yield the shifts in the uncertainties due to measuring x, with this smoothed-square- 
wave detector: 

(Ap)| = (A P ) 2 + r,(p)\ (A*)| = (Ax) 2 + a 2 x 

As you can see, for any value of the dimensionless parameter, ab, the smoothed-square-wave 
detector satisfies rj{p) 2 > 1/4<t 2 . 

4.2. Arthurs-Kelly joint measurements 

The notion of a POVM generalizes to measure spaces, X, more complicated than X C R. 
The most obvious generalization is to joint measurements, as discussed by Arthurs and Kelly 
[9]. Here, X = Spec(x) x Spec(p) = R 2 . Define 

L xoM : ' ^ J dy h(x - x ,y - x ) e ik °^ t/>(y) (38) 

where h(u,v) is normalized 

1 = J dudv \h(u,v)\ 2 
so that F XQ) k = L XQ ko L XO)ko obeys the completeness relation 

Ztx 

and therefore defines a POVM on X = Spec(x) x Spec(p) via 

^ M- Js = / ^ ,fc (39) 

Since we'll have recourse to them, let us denote the moments of this probability distribution 
by 

kl\\ f dx dk 0kij 



W= / -^xXTr(F xoM p) 



dudvdx h(u, v) (x — v) 



; d_\ l 

ds 



(40) 

(h(u, v + s)K(x + s, x)) 

=0 
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Arthurs and Kelly's choice corresponds to 

h(u,v) = g(u)J(v) (41) 

where f(v) is a Gaussian 

f ^WW^ v * (42) 

and g(u) is an (a-priori independent) acceptance function, normalized as in fl20l) . The re- 
sulting F XQ>ko , 

F XoM : V(z) H- / *//(x - x )7(x - x )e* fco(a; - y V(y) (43) 



is independent of the choice However, the density matrix, after the joint measurement, 

f dx dk t 
P - J 2tt °' k ° P x °' k ° 

does depend on that choice. If K(x, y) is the integral kernel representing p, then the integral 
kernel representing p is 



K(x, y) = I dudv h(x — v , x + u — v)h(y — v,y + u — v)K(x + u,y + u) 

(44) 

dudv g(x — v)g(y — v)f(x + u — v)f(y + u — v ) K(x + u,y + u) 



A peculiar feature of ( )4~Tj) is that the F XQjko , in ( )43|) . are projection operator^]. Moreover, 
when g = f, the L XOi fc are also projection operators. 

The measured uncertainties in x and p, for this joint measurement, are 

(Assured = «* 2 )) " W (^measured = «P 2 )) " fpf 

= (Ax)J + a 2 = (Ap) 2 + al ^ 

where 

o^ = / cW|/( M )| 2 

(46) 

a 2 = - / duj(u)f"(u)= [ du\f'(u)\ 2 



and, as usual, we've assumed that f(u) = f(—u). The Arthurs-Kelly choice (1421) of a 
Gaussian, yields a 2 = 6 2 /2 and a 2 = 1/26 2 and hence 

b 2 1 

(Az)L»ured = ( AX )1 + ^ » (Ap)Lasurcd = ( A ?)p + ^ (47) 

So, for the joint measurement, we have 

1 1 b 2 

( Aa; )mcasurcd( A P)measurcd > g + ( Ax )J &2 + ( A ^ ^ (48) 

> 1 



5 But note that the F XoM are not orthogonal: F Xo ^ F x > oik ' o ^ 0, so is still a POVM, not a PVM. 
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where the last inequality is saturated for a pure state, p, which is a Gaussian wave packet 
of very carefully chosen width: 



(Ax)J = -, (Ap)J = ^ 

Comparing with ( 13 7p . we see that this is akin to the limit a p — > for the successive mea- 
surements (so that the measured uncertainties (J25]) and (|25|) approach those of (T4T|) ). 

It is illuminating to compare the state after the successive measurements with that after 
the joint measurement. If we first measure x with the Gaussian detector (1291) and then 
measure p with the Gaussian detector (|32|) . the final state, p& na x, is given by the integral 
kernel 



K &naX (x,y) = e 



-(x-y) 2 /8a : 



2 O r , 



VW2 



du e~ 2a p u2 K(x + u,y + u) 



(49) 



which suppresses the off-diagonal elements, as before, and smears the near-diagonal elements 
against a Gaussian. This yields the net disturbance (as a result of both measurements) 



4(7?' 



rj(xy 



Aa 2 



(50) 



For the joint measurement, the final state (I44p depends on the choice of g(u), which we have 
heretofore left arbitrary. If we choose a Gaussian 



9W 



then 



becomes 



-u 2 /2a? 



K(x, y) = exp 



a 2 + b 2 



Aa?b 2 

so that the disturbances 



{x-yf 



du 

TP + 1> 2 ) 



TV 



exp 



u 



+ b 2 



a 2 + b 2 
2a 2 b 2 ' 



rj(x) 



a 2 + b 2 



K(x + u,y + u) (51) 



(52) 



depend on the additional parameter, a. 

One reasonable way to compare successive versus joint measurement, is to fix the dis- 
turbances to x and p (i.e., set the final states, (149j) and ( 15 ip . equal), and then compare the 
measured uncertainties. This is possible when the detector resolutions for the successive 
measurements satisfy a x a p < \: 



l - (l +yfl - 16a 2 a 2 ) 
~ 16^) 



Aa 2 p 



(53) 
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V 



(Apg, int - (Ap)L C essi vc = : „V V fX-r^ k 11 < (54b) 



which yields the difference in measured uncertainties 

( Ax)J oint - ( Ax)L CC ssive = (l -y/l - 16^cr p 2 ) 2 > (54a) 

Aay p (1 +v /l - 16*** *) 

The RHS of (I54bl) ranges between — 32^2 (f° r small ct x ct p ) and (for a x a p = -g). So, in 

making a joint measurement, we trade a relatively modest deterioration in (Aa;)^ ieasured for 
a potentially large improvement in (Ap) measured . 

Finally, note that Neumark's Theorem ensures that the POVM (|39[) for the joint mea- 
surement can be realized as a PVM in a larger Hilbert space, %®T-L' (in our case, we could 
choose H' = £ 2 (R 2 )). In this larger Hilbert space, the observables corresponding to this 
PVM commute, as can be seen from (1401) . Hence there exist, in this larger Hilbert space, 
states with arbitrary low values for the product of the corresponding quantum-mechanical 
uncertainties. For any (H,A,B), it is always possible to find such an augmented Hilbert 
space, and pair of commuting observables thereon. Di Lorenzo [I] would like to interpret 
this freedom as the ability to make arbitrary-precision measurements of the original (non- 
commuting) observables, A, B. But the requisite detailed knowledge of the initial state of 
the system and the ability to "prepare" the measuring apparatus (described by H') in a 
suitably corresponding initial state, belies the interpretation of this as a "measurement." 



5. Relation with Ozawa 



5.1. The disturbance rj(B) and Ozawa's uncertainty relation 

Ozawa [TlfTO] provides a more baroque definition of the disturbance, f]{B), than the one 
presented in (TTTT) above. Introduce an auxiliary Hilbert space, W and a density matrix, \ on 
it. Let U be a unitary operator on % ® H', and let the above data (H', X-, U) be constrained 
by demanding that the density matrix p, after measuring A, be given by the partial trace, 

p = Tr w {Up® X U ] ) 

Then Ozawa's definition of rj(B) is 

V (B) 2 = T rwsW (V B <g> 1 U - B <g> l) 2 p <g> X ) (55) 

It seems rather uneconomical to introduce all this auxiliary data (H', x, U) in order to define 
a quantity which should be intrinsic to (H, p, A, B). Indeed, there isn't even a uniform 
prescription for making one choice of {%' x, U), given (%, p, A, B). 

If the number of distinct eigenvalues of A is small enough, and if the Hilbert space is 
finite-dimensional, then there is a uniform prescription that we can follow. Assume that A 
has N distinct eigenvalues, A«, with multiplicities rij. Let 

N 1 

i=i 
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which has the property that 



Let 



Tr w P t (A) x = jj, Vi 



U = exp 



N 



(A) p(A) 



j=l 



We compute 



if and only if /3 satisfies 



Pi A) pP ( i A) = Tr w (Up® X U ] ) 



2(1 - cos/3) = JV 

which has solutions for N < 4. Plugging ( 1561) and ( 1571) into (1551) . we find 

77(B) 2 = Tr [((5 -Bf + B 2 - B 2 )p 



(57) 



(58 



Even though we only computed (155]) for iV = 2, 3, 4, the result is iV-independent and perfectly 
intrinsic. We take ( 158|) as a uniform definition, for all pairs of observables A, B, on any Hilbert 
spaced 

Having defined this peculiar observable, we might well ask, what does r)(B) meanl There 
is one case where the answer is clear: 

If B = B, then it follows that (B)~ = (B) . (Note that the converse statement is not 
true.) Hence, in this case, 



V (Bf = (AB)l-(AB) 2 p 

So r)(B) adds in quadrature to the quantum-mechanical uncertainty of B in the initial 
state, (AB) , to give the quantum-mechanical uncertainty in the final state, (AB)~. 



6 For any given observable, A, we can always find, in a ad-hoc way, a triple, {%' , x, U), which realizes the 
corresponding p, and then check that (|58j) is satisfied for any choice of B. When measuring x, as in ^4.1) we 
can take 

• H' = £ 2 (R) 

• x given by the pure state with integral kernel, k(x',y') = f(x')f(y'), for some choice of acceptance 
function, f(u) = f(—u). 



• U implements a particular SL(2, R) transformation. Let 

2 = ( x' ) > 

Then U : ip(x) ^ ip(Rx). 



Setting x' = y' = xq, we reproduce and verify that (|58p holds for this case, too. 
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More generally, r](B) is unrelated to the change in the uncertainty of B (which we saw, in 
§2.2[ can even decrease). Nevertheless, rj(B) > 0. 

It is claimed in p] that rj(B) = if and only if the probability distribution for B is 
the same in the state p as it was in the initial state, p. That is obviously falsfl In the 
2-state system, we saw that, for A = J x and B = b ■ J, Tj(B) — 4= (by + b 2 z ) is completely 
independent of the state p. Take p to be an eigenstate of J x (so that p — p). The probability 
distribution for B is unaffected by the measurement of A, even though we can clearly have 
7](B) > 0. 

Ozawa [1] also proposes an uncertainty relation satisfied by f]{B). In our notation, it 
reads 

e(A) (r,(B) + (AB) p ) + (AA) pV (B) > ~ \Tr([A,B]p)\ (59) 

Here, e(A) is the "noise" in the measurement of A. For an observable, A, with a pure 
point spectrum, we can consider an ideal measuring device, with e(A) = 0. However, for 
an observable, like x, with a continuous spectrum, we need to incorporate a finite detector 
resolution and e(x) = a x , the detector resolution we introduced previously. For B = p, we 
computed in ( 13TT) that r](p) > Plugging this into flo^]) yields the uninteresting result 



This "lower bound" never comes close to being saturated. We computed a sharper lower 
bound in (j37j). 

For the 2-state system, with A = J x , B = b ■ J and p — |(1 + a ■ a), we computed 
(AA) p = l -{l-alf\ v (B) = -L(bl + blf 2 

and 

- \Tr([A, B\p)\ = ^\b y a z - b z a y \ 
Plugging these into flo^]) . we obtain 

(l-al) 1/2 (bl + bl) 1/2 >±\b y a z -b z a y \ 

which, if we recall that a - a < 1, is, indeed, satisfied, but is not the best bound one could 
obtain. In fact, 

(l-alf 2 (bl + blf 2 > \b y a z -b z a y \ 



7 The weaker statement: 

i]{B) = implies that the probability distribution for B is unaltered. 
may still be true. We have been unable to find either a counter-example or a proof of this weaker statement. 
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For the pure point spectrum case, we were able to consider the case e(A) = 0. Now we 
turn to the case of a "weak measurement" of A, which is one context where one naturally 
has e(A) > 0. 



5.2. Weak Measurements, POVMs and e(A) 

Lund and Wiseman [TT] propose a test of Ozawa's ideas, using a "weak measurement" of A. 
The notion of weak measurement introduced by Aharonov et al. [12] requires a nontrivial 
extension of standard Quantum Mechanics. However, in the case of a 2-state system (and, 
more generally, for any observable A, such that Spec(A) is a finite set), Lund and Wiseman 
can avail themselves of a much simpler definition. 

A weak measurement of A is a 1-parameter family of POVMs which interpolate between 
the PVM corresponding to A and the uniform probability distribution on Spec(A) . 

More precisely, let A have a pure point spectrum consisting of N distinct eigenvalues, Aj, 
and let Fi(9) be a family of positive-semidefinite self-adjoint operators interpolating between 

F i (0) = P$ A > and F i (0 awx ) = -il 

and satisfying 

N 
i=l 

Concretely, let f(9) be a monotonic function on [0, max ], satisfying 

/(0) = 1, f(9 max ) = 

and let F^O) = L\U with 

Li = ^ (y/N-(N-l)f{6) -y/m) 1 +vW)P! A) (60) 

Explicitly, 

F l (6) = ^{l-9(0))± + 9(d)Pi A) (61) 

where 

9(0) =( 1 ~l) + J^ N ~ ( N ~ l )WVW) (62) 
von Neumann's formula (|17[) for the density matrix after measuring this POVM is 

P = L{ P Lj i 

i 

Using fl60|) . we compute 
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p= (l-f(O))p + f(0)P 

which interpolates between p (the change resulting from an ordinary measurement of A) and 
p (no change Similarly, for any observable O, we can define 

= Y, UOL\ = (1 - f(9))0 + f(0)d 

i 

which interpolates between O and O. 

For the 2-state system, Lund and Wiseman made a specific choice, 

/(0) = l-sin(20), ma x = vr/4 

which follows from an auxiliary triple (H', x, U) as in §5.11 Their choice (with a slight change 
in conventions (a x <H> a z ), since we want A = J x whereas they chose A = 2J Z ) was 

w = u 

X = T^ 1 + cos(2#y x . + sm{26)a z ) 

U=-(l + a x )®l + -(1 - <r x ) ® a z 
Of course, there's nothing special about that choice; so long as we satisfy the constraint that 

p = Tr w (U p® X U ] ) (63) 

we can use any choice of (7i', Xi U). In fact, for N < 4, we can avail ourselves of the triple 
(H 1 , X-, U) introduced in (156]) and (jSTj) by simply allowing (3 in fl57|) to be ^-dependent 



which results in 



oiV-3 

f3(6) = 2(6 m ^-6), 9 max = —ir 



f(9) = ^(1 - 008 0(0)) 



which agrees with Lund- Wiseman for N = 2. 

The disturbance in 5, resulting from a weak measurement of A can be defined similarly 

to (El 



8 To complete the prescription of a weak measurement [13] . one needs to define a set of contextual values 
[14] . \i(6), such that 

A = ^2Xi(j3)Fi(e), V0e[O,0 maJC ) 



With our parametrization, 
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r]e(B) 2 = Tr 
A simple calculation yields 



(B-B) 



B< 



B 2 )P 



(64) 



Ve(B) 2 = f(0)v(B) 2 

with r](B) 2 given by (!58|) . We also need to introduce e(A) 2 , which represents the "noise" 
inherent in a "weak measurement" of A. 

Lund and Wiseman define e(A) 2 in terms of a quadruple (H', x, U, M): 



e(Af 



Tr 



'umv 



[U ] 1 <g> M U - A® l) 2 p ® x 



(65) 



where %' , Xi U are as before and M is some self-adjoint operator on %' . Unfortunately, aside 
from the obvious requirement that ( )63l) be satisfied, this "definition" is a little . . . under- 
specifed. Instead, we directly implement the notion that e(A) 2 represents the mean-squared 
deviation of the weakly-measured A from its true value by writing 



e{A) 2 = J> 4 - A,) 2 Tr (f^P™ p 



(66) 



i.j 



Here, Tr{Fi{9)P- pj is the probability that a "weak measurement" of A yields the value 
Xi, whereas a measurement of A would have yielded Xj. Using ( 1ST)) , this reduces to 

^) 2 = ^(l-^))£(A>-A,) 2 Tr (ifV) 

which agrees with their explicit realization of ( 165]) for the 2-state system (the only case we 
will actually need). 

Putting these together, we obtain the Ozawa/Lund- Wiseman uncertainty relation for a 
weak measurement of A = J x , followed by a measurement of B = b ■ J, starting with the 
initial state p = \ (1 + a ■ a), is: 



2 sin(0) (cos(0) - sin(0)) {b 2 y + b 2 z ) 1/2 +V2sin(0) (b 2 - (a ■ 6) 2 ) + 

1/2 1/2 1 

+ (1 - a 2 x ) (b 2 y + b 2 z )< >— \b y a z - b z a y \ 

In §5.14 we already saw that this inequality was ineffective (never close to being saturated) 
at 9 = 0. Here, we see that the situation is even worse at nonzero 6, since the "new" terms 
on the LHS are non-negative for < 9 < tt/A. 

Rozema et al. [2], in a clever experiment, measure the quantities appearing on the LHS 
and RHS of (JBTl) separately. Indeed, their results for the LHS nowhere approach the values 
for the RHS. 
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Finally, let us briefly return to the subject of §4.11 Measuring x, with finite detector 
resolution a x , is naturally associated to the POVM 



and a simple computation yields 
as expected. 
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The analogue of (1661) is 




(68) 
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